Universal Reduction of Effective Coordination Number 
in the Quasi-One-Dimensional Ising Model 
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Critical temperature of quasi-one-dimensional general-spin Ising ferromagnets is investigated by 
means of the cluster Monte Carlo method performed on infinite-length strips, LxooorLxLx 
co. We find that in the weak interchain coupling regime the critical temperature as a function 
of the interchain coupling is well-described by a chain mean-field formula with a reduced effective 
coordination number, as the quantum Heisenberg antiferromagnets recently reported by Yasuda et 
al. [Phys. Rev. Lett. 94,217201 (2005)]. It is also confirmed that the effective coordination number 
is independent of the spin size. We show that in the weak interchain coupling limit the effective 
coordination number is, irrespective of the spin size, rigorously given by the quantum critical point 
of a spin-1/2 transverse-field Ising model. 



I. INTRODUCTION 

Effects of interchain coupling in quasi-one-dimensional 
(Q1D) magnets have been extensively investigated over 
many years. Theoretically, introduction of an infinites- 
imal interchain coupling can bring a drastic change to 
the system; emergence of a thermal phase transition and 
long-range magnetic order therewith at finite tempera- 
tures. In interpreting experimental results on Q1D ma- 
terials, which mostly exhibit anomalous behavior at fi- 
nite temperatures, one often has to take the interchain 
effects into account as well. A standard and established 
approach for incorporating such a weak interchain inter- 
action is the chain mean-field (CMF) approximation, 1 - 2 
in which a problem in two or three dimensions is reduced 
to a single-chain model in an effective field by ignoring 
fluctuations of order parameter between weakly coupled 
chains. Within the CMF theory, the magnetic suscepti- 
bility of a Q1D system is approximated as 



X(T) 



Xid(T) 



l-zJ'xw(T) 



(1) 



where Xid{T) is the susceptibility of a genuinely one- 
dimensional chain, z the coordination number, i.e., the 
number of nearest neighbor sites in directions perpen- 
dicular to the chain, and J' is the interchain coupling 
constant in these directions. Since Xid{T) may not have 
any singularity at finite temperatures, the critical tem- 
perature T c of the Q1D system is given as the locus of 
the simple pole in the right-hand side (rhs) of Eq. QJ. 
That is, 



1 



J'Xu(T c ) 



(2) 



is the equation which determines the Q1D critical tem- 
perature. It has been confirmed that for various Q1D 
(and also Q2D) materials the solution of the CMF self- 
consistent equation J5J) reproduces a qualitatively accept- 
able J'-dependence of the critical temperature in the 
weak interchain coupling regime. 



Recently, Yasuda et al£ carefully reexamined the Neel 
temperature of the quantum Heisenberg antiferromag- 
nets on a anisotropic cubic lattice, and found that Eq. @ 
does not hold even in the weak interchain coupling limit. 
In Ref. 11 the rhs of Eq. Q for various J"s is evalu- 
ated accurately by the extensive quantum Monte Carlo 
simulation using the continuous-imaginary-time loop al- 
gorithm, and it is concluded that for J' — > it converges 
not to 4 but to 2.78 for the Heisenberg antiferromag- 
nets on the Q1D cubic lattice. This discrepancy with the 
CMF theory is interpreted as a reduction of the effective 
coordination number 



a j') = 



i 



J'xidW)) 



(3) 



from that of the original lattice, z. Furthermore, the 
numerical data presented in Ref. clearly demonstrate 
that the renormalized coordination number is "univer- 
sal," i.e., independent of the spin size S (including the 
classical limit S — oo). Although there have been sev- 
eral attempts^'SiS'^iilliii to develop theories beyond the 
CMF for specific models so far, our understanding of S- 
independent renormalization of the effective coordination 
number is, fairly speaking, not yet satisfactory. 

Another interesting question is whether and how such 
a renormalization of the effective coordination number is 
observed for other models than the Heisenberg antiferro- 
magnets. Among many theoretical spin models, a linear 
array of S — 1/2 Ising chains, which forms an anisotropic 
square lattice, provides a valuable example, since the ex- 
act critical temperature is known for any ratio of J' to 
the intrachain coupling J. Indeed it was pointed out 
more than three decades agoi that the CMF result for 
the critical temperature 

T ~ - (4) 

c 21n(J/J') - 2 In ln( J/J') 1 ' 

for J'/ J -C 1 is not compatible with the exact asymptotic 
behavior 



J 



21n(2J/J') - 21nln(2J/J') 



(5) 
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Here and hereafter we put the Boltzmann constant to 
unity. The former asymptotic expression is obtained from 
Eq. (J2J with z = 2, combined with the exact susceptibil- 
ity of an 5 = 1/2 single Ising chain 



(a) 



(b) 



J_ J/2T 

AT 



(6) 



whereas the latter is obtained by expanding Onsager's 
exact solution for the anisotropic Ising models 



sinh — — sinh — — 
2T r 2T r 



(7) 



for small J 1 / J. [Note that here the absolute spin length 
is not 1 but 1/2. See the definition of our Hamiltonian (JSJ 
below.] Interestingly, if the effective coordination num- 
ber [Eq. is evaluated for this case by using the ex- 
act critical temperature and the exact one-dimensional 
susceptibility, Eqs. and J7J) respectively, one immedi- 
ately finds that it converges to unity, instead of 2, in the 
J'/ J = limitiiS This example thus clearly demonstrates 
the breakdown of the CMF approach for Q1D magnets, 
to which we will address ourselves in the present work. 

In the present paper, we report our detailed analysis 
on the Q1D general-spin Ising models, and discuss the 
origin of the 5-independent renormalization in the effec- 
tive coordination number in the weak interchain coupling 
regime. The organization of the paper is as follows: In 
Sec. [n] we define our models and summarize some ba- 
sic properties in the genuinely one-dimensional case. In 
Sec. 11111 we describe our Monte Carlo method and the 
detail of finite-size scaling analysis. The main results 
of our Monte Carlo simulation are presented in Sec. IIVI 
where we unambiguously show that the critical temper- 
ature is described quite well by a CMF formula with a 
reduced effective coordination number. The effective co- 
ordination number is found independent of 5 in the weak 
interchain coupling limit, but does depend on the dimen- 
sionality of the lattice as well as on the fine structure 
of the lattice in directions perpendicular to the chains. 
In Sec. [V] we discuss the origin of the "universality" (or 
"nonuniversality" ) by considering a mapping to a quan- 
tum Ising model in the weak interchain coupling limit. 
We show that the original spin-5 Ising model is mapped 
onto the spin- 1/2 transverse-field Ising model, irrespec- 
tive of the spin size, and furthermore the effective coordi- 
nation number in this limit is rigorously independent of 5 
and given by the quantum critical point of the mapped 
5 = 1/2 transverse-field Ising model. Our results and 
conclusions are summarized in Sec. I VII 



II. QUASI-ONE-DIMENSIONAL ISING MODEL 

We consider the following classical spin-5 Ising model: 

i.r i.r,8 
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FIG. 1: Lattice structures in the transverse directions 
perpendicular to the chains: (a) square lattice (SQ), (b) 
square lattice with anisotropic next-nearest-neighbor interac- 
tion (SQ+ANNNI), and (c) simple cubic lattice (SC). All the 
bonds shown in the figure have the exchange coupling of the 
same strength, J'. 



where i and r are the position of lattice sites along chains 
and perpendicular directions, respectively, 8 is summed 
over the z/2 nearest- neighbor vectors in the transverse 
directions, and Sf r (— S, 5 — 1, . . . , — 5) is the spin vari- 
able at site (i, r). Both of the longitudinal and transverse 
coupling constants, J and J', are assumed to be ferro- 
magnetic (J > J' > 0). 

In the present study, three different Q1D lat- 
tices are considered: square lattice (SQ), square lat- 
tice with anisotropic next-nearest-neighbor interaction 
(SQ+ANNNI), and simple cubic lattice (SC). In Fig. [TJ 
the lattice structure in a slice perpendicular to the chains 
are drawn. The former two (SQ and SQ+ANNNI) are 
two-dimensional lattices, while the last one (SC) is three 
dimensions. On the other hand, the latter two shares the 
same coordination number, z = 4, which is larger than 
that of the first lattice (z = 2). By considering these 
three lattices, one can evaluate separately the effects of 
dimensionality from the coordination number. 

The linear dimension of the lattice in the directions 
perpendicular to the chains is specified by L. We assume 
periodic boundary conditions in these directions. For the 
SQ and SQ+ANNNI lattices, the number of spins in each 
slice, denoted by N s hereafter, is equal to L, while it is 
L 2 for the SC lattice. In the chain direction, on the 
other hand, the system size is assumed to be infinite as 
explained in the next section. 

Before starting a detailed discussion about the criti- 
cal temperature of the Q1D system, let us briefly review 
thermal properties of a genuinely one-dimensional single 
chain, i.e., Eq. JSJ with J' = 0, since its susceptibility 
Xid(T) appears as a key quantity in the CMF approx- 
imation, and thus is expected to dominate asymptotic 
behavior of the critical temperature T C (J') for J' « J. 

The partition function of a periodic spin-5 Ising chain 
of length M is represented asi£ 



trT 



M 



(9) 



Here the transfer matrix T is a (25 + 1) x (25 + 1) real 
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symmetric matrix with matrix elements 

T .. = e J(S+l-i)(S+l-j)/T_ (1Q) 

The matrix T is diagonalized as 



T = UDU 



(11) 



by a diagonal matrix of eigenvalues, = Aj<5j ;J -, and a 
matrix of the corresponding right (column) eigenvectors, 
U = (u%, U2, . • • , "25+1) with Tui — XiUi. We assume all 
the eigenvectors are normalized, i.e., U^ 1 = t U. With- 
out loss of generality we further assume Ai and A2 denote 
the largest and the second largest eigenvalues, respec- 
tively. 

At finite temperatures, the transfer matrix T is pos- 
itive definite and its largest eigenvalue Ai is nondegen- 
erate. Thus the partition function is evaluated as 
(Ai) M in the thermodynamic limit. Similarly, the ther- 
mal average of correlation function between the spins at 
site and r is calculated in terms of the transfer matrix 



(S*S?) = -tvST r ST 



M-r 



= \tv{U- l SU)D r {U- l SU)D 
Af r tr {U- 1 SU)D r (U~ 1 SU) 



M-r 



(12) 



2S+1 



= ^ \ tu l Su i 



where we introduce Si 



Ai 



(5+1 



i)§i,j- In the last 
= 1 in the summa- 



expression, we explicitly dropped i 
tion, since (*iti5ui) constantly vanishes due to the spin 
inversion symmetry in the Hamiltonian JSJ) . The suscep- 
tibility is then obtained by summing up the correlation 
function over the whole lattice 



Xid(T) = - 



E 



(S Q S r ) — 



1 25+1 x _i_ \ 

It C< |2 A l + A l 

^ I mS Ui \ -— -- 

i=2 1 1 



(13) 



For 5 = 1/2, 1, and 3/2, all the eigenvalues and the 
eigenvectors of the transfer matrix can be obtained ex- 
plicitly in a closed form, and so is the susceptibility^ For 
larger 5, it is not able to solve the eigenvalue problem 
analytically, but one can still diagonalize T numerically 
very easily. 

At very low temperatures, T <C J, the largest two 
eigenvalues Ai and A2 nearly degenerate with each other, 
which reflects two degenerating fully polarized ground 
states, (5, S, . . . , S) and its spin-inverted counterpart, 
with energy —S 2 JM. The corresponding eigenvectors 
are symmetric or antisymmetric, respectively, under the 
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FIG. 2: Scaling plot of susceptibility of a single chain with 
5 = 1/2 (solid line), 1 (diamonds), 3/2 (upward triangles), 2 
(downward triangles), 4 (circles), and 8 (squares). In the 
inset, corrections to the low-temperature asymptotic form 
[Eq. 1241 1 are presented. 



spin inversion, and their elements are given as follows: 

(14) 



U2.1 



1 = tti,asr+i = A= + 0(e J / T ), 



0{e~ J ' T ) fori = 2,..., 25, 
1 



'"1,25+1 



x/2 



0(e 



u 2li = 0(e~ J/T ) for i = 2, . 



-J/T\ 



,25. 



(15) 
(16) 
(17) 



The absolute value of the subdominant eigenvalues is ex- 
ponentially small in comparison with the largest two. 
The correlation function and the susceptibility at low 
temperatures are thus written as 



(5 z 5*) ~ 5 2 (^) r and 
5 2 Ai + A 2 



Xid(T) 



T Ai - A 2 



(18) 
(19) 



respectively. By defining the correlation length £id (T) as 
Zid( T ) = -log(A 2 /Ai), we obtain 

5 2 1 + e-VfidCO 25 2 
Xi d (T).- T ±- T7? - M .— a d (T) (20) 

for general spin 5. For 5 = 1 and 3/2, the low- 
temperature susceptibility is explicitly given by«i^ 



X \ d (T) ~ ±e^ T and 



X% 2 (T) 



° 9J/2T 

4T 



(21) 
(22) 



respectively. [The exact expression for 5 = 1/2 is already 
presented in Eq. ©.] For larger 5, analytic expressions 
are not known, but low-temperature asymptotic forms 



4 



for the correlation length and the susceptibility are con- 
jectured as follows: 



&(T)~- 251og(tanh^ 



J_ 2S 2 J/T _ }_ flrp C 2 

45 



T 

f(T/S*J) 



(23) 



and 



xfdCO 



(24) 



where f(x) — e 2 l x /4 and g(x) — e 2 ' x /2x are scaling func- 
tions. These expressions are confirmed explicitly, e.g., 
Eqs. ©, iJHJ, and $22$, for 5 < 3/2. We emphasize 
that these scaling forms are not very trivial. Indeed, a 
naive dimensional analysis yields £f d (T) — f(T/S 2 J) and 
xf d (T) ~ J~ 1 g(T/S 2 J) with some scaling functions f(x) 
and g{x), which differs from the correct scaling forms 
presented above by a factor 1/5. In Fig. [3 we present 
a scaling plot of the numerically exact data for the sus- 
ceptibility, given by Eq. for 5= 1/2, 1,.. .,8, which 
confirms the asymptotic scaling formula (|24(l . 



III. NUMERICAL METHOD 

A. Cluster Monte Carlo on infinite-length strip 

For the 5=1/2 SQ lattice, the critical temperature 
is solved exactly for any combination of J' and / J as 
we have already seen in Sec. [I] However, no exact so- 
lutions are available for 5 > 1. For the other lattices 
(SQ+ANNNI and SC) there are no solutions even for 
5 = 1/2. 

In the present study we employ the Monte Carlo 
method combined with the finite-size scaling analy- 
sis for evaluating unbiased critical temperatures of the 
Q1D Ising models. Since the present model has only 
the ferromagnetic couplings and has no frustration, the 
Swendsen-Wang cluster algorithm^ or the Wolff cluster 
algorithm^ is one of the best algorithms. To optimize 
further for the Q1D general-spin Ising models, we extend 
the cluster Monte Carlo algorithm in the following two 
points. 

The first optimization is for strong anisotropy. Since 
our model is strongly spatially anisotropic, one has to 
fine-tune the aspect ratio of the lattice to simulate an 
effectively large system with a reasonable costi&ii In the 
present study, instead of anisotropic finite lattices, we 
simulate infinite length strips N s x oo directly, follow- 
ing the infinite-lattice algorithm proposed by Evertz and 
von der Linden?i& In the original proposal of the infinite- 
lattice cluster algorithm, the lattice is assumed to be in- 
finite in all the spatial directions and only one cluster 
including the origin is build at each Monte Carlo step. 
In our simulation, however, the lattice is infinite only in 
one direction, the chain direction, which has a stronger 



coupling than the transverse directions, and we build up 
completely all clusters including the N s sites on the slice 
in the center of the lattice at each Monte Carlo step. 
This modification avoids the divergence of cluster size 
and enables us the simulation at and below the critical 
temperature of the bulk and thus the finite-size scaling 
for accurate estimation the critical temperature. 

Another modification to the original cluster algorithm 
is the extension for general spin 5 > l/2m&2& To this 
end, first we decompose a spin into 25 subspins, that is, 
we simulate the following 5=1/2 model 



2S 

Wsubspin = — J y ] y ] ^i,r,a^i+l,r,a' 
i,r a,a' = l 
2S 

- T' \^ \^ q z S! z 

° / i / j u i,r,a u i,r+8,a'i 



(25) 



instead of the original spin-5 model. Here 5j irjQ , = ±1/2 
is the ath subspin at site (i, r). The volume of the phase 
space of the mapped system l|25l) is enlarged than the 
original one [from (25 + 1) N to 2 2SN , where N is the 
number of original spins] . We need an extra procedure 
which symmetrize each spin to recover the original phase 
spac o 1 ?' 2 ? as described below. 

A Monte Carlo step is as follows: (i) set i i— > 0; (ii) 
every possible subspin pair, Si. r , a and Si ir _|_ s,a'i in the 
ith slice is connected with probability (1 — e~ 2J / T ) if 
Si, r ,a — Si tr _|_ s,a' , or left disconnected otherwise; (iii) 
for each group of subspins (i,^), the subspins are parti- 
tioned into two according to their spin direction, a ran- 
dom permutation is generated for each part, and the 
subspins belonging to the same permutation cycle are 
connected with probability one; (iv) every possible sub- 
spin pair, 5i i7 - jQ and 5i4.i jI%Q ', on the neighboring slices, 
i and i + 1, is connected with probability (1 — e~ 2J / T ) 
if Si, rt0t = 5j + i irjQ /, or left disconnected otherwise; (v) 
increase i by 1 and repeat (ii)-(iv) until no clusters in the 
current slice include any subspins on the central slice at 
i = 0; (vi) repeat (ii)-(v) in the opposite chain direction, 
i = —1,-2,...; (vii) generate a new spin configuration 
by flipping all subspins on each cluster including the sub- 
spins on the central slice randomly with probability 1/2, 
and discard all the clusters. 

The present method satisfies detailed balance condi- 
tions and ergodicity for any finite region including the 
central slice. Typically, we spend 1.0 x 10 6 Monte Carlo 
steps for measurement after discarding 3.3 x 10 4 steps for 
thermalization for each parameter set (J' / J,T/ J, L). 



B. Susceptibility and correlation length 

Correlation functions {S* r S z , r ,) can be measured effi- 
ciently in terms of a neat technique, so called improved 



5 



estimator^ i.e., 

(Sl r Sf, y ) = i]T(C(z,r,a;z',r>')}MC, (26) 



where (• • • )mc denotes a Monte Carlo average and 
C(i, r, a; i', r', a') takes 1 if the subspins at (i, r, a) and 
(i',r',a') belong to the same cluster and takes other- 
wise at each Monte Carlo step. We can further rewrite 
Eq. (|26|l by introducing a function Ce(i, r, a), which takes 
1 if the subspin S' ra is on the £th cluster and takes 
otherwise, as 



4 ^ 



^2Ct(i,r,a)C e (i',r', 



a 



MC 



(27) 

Here one must recall that in the present algorithm clus- 
ters will be built only for a part of the lattice, and there is 
no way to judge whether two subspins belong to the same 
cluster when neither is on the built clusters. Therefore 
the above improved estimator for the correlation func- 
tion is valid only if at least one of (i,r,a) and (i', r',a') 
is guaranteed to belong always to built clusters. Since 
only the subspins on the central slice satisfies this condi- 
tion, i or i' must be zero in Eqs. and (|2T|) . 

The magnetic susceptibility \ is obtained by integrat- 
ing the correlation function over the whole lattice. Fix- 
ing i = in Eq. (|26[1 and taking a sum over all the other 
indices, we obtain the following expression for the sus- 
ceptibility: 



TN, 



OO 



Z— — 00 r,r' 



1 OO 

W £ £ (£<7H0,r,a)Q(i,r->'; 



4T/V S 



%— — oo r .r'a.a' 



ATN f 



MC 



MC 



(28) 



where nf ^ = J2 r a Q(®> r > a ) denotes the number of 
subspins on the central slice in the £th cluster, and 
ri£ = Ci(i, r, a) is the number of subspins in the 

tth cluster. 

The correlation length along the strip is estimated as 



X k \l/(k-k') 



(29) 



with nonnegative integers, k ^ k', and the moments of 
the correlation function 



oo 



i— — oo r.r 



4fc!/V s \^ e Vmc' 



(30) 
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FIG. 3: Scaling plots of the correlation length for the S = 1/2 
SQ lattices with J'/J = 1 (left) and J'/J = 0.1 (right). 
System sizes are L = 22, 24, ... , 30 and L = 14, 16, ... , 24, 
respectively. The error bar of each point is much smaller 
than the symbol size. The critical temperature and the crit- 
ical amplitude are estimated as T C /S 2 J = 2.2695(3) and 
L/£\t=t c = 0.789(3) for J'/J = 1, and T C /S 2 J = 0.9058(2) 
and Lji \ t =t c = 0.176(1) for J'/J = 0.1, respectively 



where m k = \i\ k Ce(i, r, a). Note that the moment 

Xk converges for any integer k > 0, since the correlation 
function decays exponentially on a finite- width strip even 
below the critical temperature. In the following analysis 
we use the simplest estimator, £1,0 = Xi/Xq, since its 
statistical error is the smallest among all the combina- 
tions (k, k') we examined. 



C. Finite-size scaling analysis 

Near the critical point T = T c the correlation length 
and the susceptibility of a finite system obey the following 
finite-size scaling formulas: 

L/ttL,T) = ttL l/v {T-T c )/T c ) and (31) 
L-^ X (L, T) = x(L 1/u {T - T c )/T c ), (32) 

where v and 7 are the critical exponents of the correla- 
tion length and of the susceptibility respectively and £ (x) 
and x(x) are scaling functions. In the following analysis 
we assume the critical exponents of the Ising universal- 
ity, i.e., v = 1 and j/u = 7/4 for SQ and SQ+ANNNI, 
and v = 0.6289 and j/u = 1.9828 for SC. 2 i The critical 
temperature is then determined as the point where all 
the data points collapse best on a single curve. 

As a test of the numerical method explained in this 
section, we perform a Monte Carlo simulation and a 
finite-size scaling analysis for S — 1/2 SQ lattices with 
J'/J = 1 and 0.1. In Fig. [31 we presents scaling plots 
of the correlation length. All the data points collapse 
excellently on a single curve in both cases. The criti- 
cal temperature is estimated as T C /S 2 J = 2.2695(3) and 
0.9058(2) for J'/J = 1 and 0.1, respectively which coin- 
cide with the exact values T C /S 2 J « 2.26919 and 0.90588. 
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Here we comment that the present lattice geometry 
L x oo is the same as the one used for the phenomeno- 
logical renormalization calculation of the transfer ma- 
trix. As a result, not only the critical exponents but also 
the critical amplitude of the correlation length along the 
strip is given analytically in the isotropic case (J' = J) 
due to the conformal invariance of the model at the 
critical points We confirm that the critical amplitude 
L/£\ t =t c = 0.789(3) for J' /J = 1 (Fig. 01 agrees with 
7r/4 « 0.785, the prediction from the conformal field the- 
ory. In anisotropic cases J' ^ J, a nonuniversal pref ac- 
tor, which represents effective aspect ratio, is introduced 
and thus the overall critical amplitude becomes nonuni- 
versal. However, this gives us a chance to estimate the 
effective aspect ratio at the critical point directly. For ex- 
ample, we can estimate the effective aspect ratio immedi- 
ately as 0.789/0.176 w 4.48 for J'/ J = 0.1. This method 
works also for the three-dimensional case (but does not 
for the SQ+ANNNI lattice due to the anisotropy intrin- 
sic to the lattice). On the other hand, it is practically a 
hard task if one works in the standard square or cubic 
geometry, L x V or L x V x L", for Monte Carlo simu- 
lations. This is another advantage of the present Monte 
Carlo method for analyzing spatial anisotropy precisely, 
though we will not utilize this quantity in the remainder 
of this paper. 



IV. RESULTS OF MONTE CARLO 
SIMULATION 

By using the Monte Carlo method and the finite-size- 
scaling technique explained in the last section, we esti- 
mate the critical temperature of the SQ, SQ+ANNNI, 
and SC lattices very carefully for 5 — 1/2, 1, and 3/2 
and J' I J = 0.001, 0.002, 0.005, . . . , 1. The error bar of 
the estimates is typically less than 0.01%. Such a high 
accuracy for the critical temperature is essential for a 
reliable estimation of the effective coordination number 
especially for small J'/ </, where any small error in T c 
is amplified vastly due to the exponential growth of the 
single-chain susceptibility at low temperatures. The re- 
sults are summarized in Tables [H] and IIIII together 
with the maximum system size for each value of J'/ J. 

The present estimates of the critical temperature for 
the 5=1 and 3/2 isotropic SQ lattices (J 7 / J — 1) can 
be compared with the previous high-temperature series 
expansion study,^ T C /S 2 J = 1.69356(1) and 1.46144(1), 
respectively. The result for the 5 = 1/2 isotropic SC 
lattice also agrees with the previous Monte Carlo study, 21 
T c /5 2 J = 4.511424(53) within the error bar. 

For all the lattices and all the spin sizes we examined, 
the critical temperature, normalized by J, becomes lower 
with decreasing J'/ J, as expected. The decrease is found 
extremely slow. Even at J'/ J = 0.001 the critical tem- 
perature is as high as 10%-20% of the isotropic case. This 
slow decrease is explained, within the CMF theory, as a 
reflection of exponential divergence of the single-chain 
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FIG. 4: J'-dependence of the critical temperature for the SQ 
lattice by the CMF approximation (dashed lines) compared 
with the exact [S — 1/2 (squares)] and the Monte Carlo [S = 
1 (diamonds) and 3/2 (circles)] results. Solid lines denote the 
CMF results with a renormalized coordination number f = 1 
(see the text for details). The error bar of the data for S — 1 
and 3/2 is much smaller than the symbol size. 



susceptibility [Eq. JHJI]. 

In Fig. 21 we plot the critical temperature of the SQ 
lattice as a function of J'/ J together with the estimates 
by the CMF approximation, i.e., the numerical solutions 
of Eq. (J2J. One sees that the CMF theory constantly 
overestimates the critical temperature. This tendency is 
naturally understood as a systematic error subsistent in 
the approximation, introduced by ignoring fluctuations 
of order parameter between chains. Although the devia- 
tion of the CMF results from the correct T C (J') becomes 
smaller with decreasing J'/ J, its convergence is logarith- 
mically slow for all the spin sizes, as demonstrated for 
5 = 1/2 in Sec. |U Similar J'-dependence of the critical 
temperature is observed for the other lattices. 

In Sec.[IJ we see that the effective coordination number 
for the 5=1/2 SQ lattice converges to unity instead of 
2, which indicates that the critical temperature in the 
weak interchain coupling regime is well-described by a 
CMF approximation with the renormalized coordination 
number, £ = 1. In Fig. we also plot the results of 
the CMF approximation with ( = 1 by solid lines, which 
amazingly agrees with the true critical temperature in a 
very wide range of J' / J; even at J'/ J = 1, deviation from 
the correct value is only 3%, 2%, and 4% for 5 = 1/2, 1, 
and 3/2, respectively, which is much smaller than those 
of the conventional CMF theory, 55%, 43%, and 39%. 
For J' / J < 0.1, deviation cannot be recognized anymore 
by the present scale of the vertical axis in Fig. 0] 

To investigate the coordination number reduction in 
more detail, we next plot the effective coordination num- 
ber [Eq. |0|] as a function of J'/ J in Fig. |3J Precise 
values are also listed in Tables HJ [HI and IIIII One sees 
in Fig. [3] that the effective coordination number clearly 
converges to a finite value very rapidly for all the lat- 
tices and all the spin sizes. The 5-independence of the 
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TABLE I: Critical temperature and effective coordination number of the SQ lattice for 5 = 1/2, 1, and 3/2. The critical 
temperature for S = 1/2 is obtained from Eq. 0. The effective coordination number for 5 = 1 and 3/2 in the J'/J = limit, 
presented in the bottom row, is estimated by taking an average for J' / J = 0.001, 0.002 and 0.005. 







5 = 


1/2 


5 = 


= 1 


5 = 


3/2 


J j J ^max 


T C /S 2 J 




T C /S 2 J 




T C /S 2 J 




l 


72 


2.26919 


0.939927 


1.6934(2) 


1.0492(2) 


1.4611(1) 


1.0907(2) 


0.5 


48 


1.64102 


0.970162 


1.2352(1) 


1.0430(2) 


1.0686(1) 


1.0731(3) 


0.2 


44 


1.14156 


0.989893 


0.88657(5) 


1.0230(2) 


0.77513(5) 


1.0393(2) 


0.1 


40 


0.905883 


0.995958 


0.72490(7) 


1.0111(4) 


0.64059(6) 


1.0191(4) 


0.05 


36 


0.741313 


0.998486 


0.61061(3) 


1.0045(2) 


0.54589(3) 


1.0082(3) 


0.02 


36 


0.590681 


0.999618 


0.50250(2) 


1.0010(2) 


0.45587(2) 


1.0021(2) 


0.01 


36 


0.508926 


0.999871 


0.44146(2) 


1.0002(3) 


0.40448(2) 


1.0004(3) 


0.005 


32 


0.445462 


0.999958 


0.39256(2) 


0.9998(3) 


0.36285(2) 


0.9999(4) 


0.002 


32 


0.380981 


0.999991 


0.34134(1) 


0.9999(2) 


0.31863(2) 


0.9997(5) 


0.001 


32 


0.342657 


0.999997 


0.31008(2) 


0.9997(5) 


0.29131(2) 


1.0000(6) 








1 




0.99989(6) 




0.99986(9) 



TABLE II: Critical temperature and effective coordination number of the SQ+ANNNI lattice for 5 = 1/2, 1, and 3/2. The 
effective coordination number in the J'/J = limit, presented in the bottom row, is estimated by taking an average for 
J'/J = 0.001, 0.002 and 0.005. 







5 = 


1/2 


5 = 


= 1 


8 = 


3/2 


J J J -£/max 


T C /5 2 J 


C 


T C /S 2 J 


c 


T C /S 2 J 


C 


l 


72 


3.8487(9) 


2.2889(8) 


2.8064(4) 


2.5325(6) 


2.4032(2) 


2.6157(3) 


0.5 


48 


2.5799(3) 


2.3766(5) 


1.8790(3) 


2.5656(8) 


1.6084(2) 


2.6223(6) 


0.2 


44 


1.6451(1) 


2.4388(3) 


1.2240(1) 


2.5495(5) 


1.0548(1) 


2.5947(6) 


0.1 


40 


1.2378(1) 


2.4600(5) 


0.9478(2) 


2.522(2) 


0.8244(1) 


2.550(1) 


0.05 


36 


0.97023(6) 


2.4698(5) 


0.76744(5) 


2.4990(6) 


0.67517(7) 


2.513(1) 


0.02 


36 


0.73958(5) 


2.4747(6) 


0.60904(4) 


2.4831(7) 


0.54444(4) 


2.4889(8) 


0.01 


36 


0.62076(3) 


2.4757(5) 


0.52440(4) 


2.4786(9) 


0.47414(3) 


2.4810(8) 


0.005 


32 


0.53183(3) 


2.4751(7) 


0.45878(2) 


2.4768(6) 


0.41909(2) 


2.4768(7) 


0.002 


32 


0.44467(2) 


2.4756(6) 


0.39197(1) 


2.4764(4) 


0.36232(2) 


2.4755(9) 


0.001 


32 


0.39441(2) 


2.4757(8) 


0.35215(1) 


2.4758(5) 


0.32803(2) 


2.476(1) 








2.4755(2) 




2.4763(3) 




2.4763(4) 



TABLE III: Critical temperature and effective coordination number of the SC lattice for 5 = 1/2, 1, and 3/2. The effective 
coordination number in the J'/J = limit, presented in the bottom row, is estimated by taking an average for J'/J = 0.001, 
0.002 and 0.005. 

5 = 1/2 5=1 5 = 3/2 



J j J ^raax 


T C /S 2 J 


C 


T C /S 2 J 


c 


T C /S 2 J 


c 


l 


22 


4.5117(3) 


2.8962(3) 


3.1966(4) 


3.0814(6) 


2.7138(3) 


3.1435(5) 


0.5 


20 


2.9297(3) 


2.9606(5) 


2.0875(1) 


3.1077(3) 


1.7746(1) 


3.1583(3) 


0.2 


20 


1.8139(1) 


3.0111(3) 


1.32750(6) 


3.1015(3) 


1.1378(2) 


3.137(1) 


0.1 


20 


1.34307(7) 


3.0296(4) 


1.01438(5) 


3.0836(4) 


0.87813(2) 


3.1060(2) 


0.05 


18 


1.03984(7) 


3.0387(6) 


0.81302(2) 


3.0654(3) 


0.71228(1) 


3.0787(2) 


0.02 


18 


0.78297(2) 


3.0434(3) 


0.63900(1) 


3.0514(2) 


0.56905(2) 


3.0563(5) 


0.01 


18 


0.65251(1) 


3.0440(2) 


0.54724(1) 


3.0469(3) 


0.49310(1) 


3.0486(3) 


0.005 


14 


0.55589(3) 


3.0444(8) 


0.47673(2) 


3.0455(7) 


0.43420(2) 


3.0459(8) 


0.002 


12 


0.46200(1) 


3.0447(4) 


0.40544(1) 


3.0445(4) 


0.37385(2) 


3.044(1) 


0.001 


10 


0.40832(2) 


3.0463(8) 


0.36325(3) 


3.046(2) 


0.33763(3) 


3.045(2) 








3.0449(4) 




3.0449(3) 




3.0453(6) 
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limiting value is confirmed to an accuracy of 0.01% (see 
the bottom row of Tables [I] |H] and . We conclude 
£ in the weak interchain coupling limit as 1, 2.4755(2), 
and 3.0449(4) for the SQ, SQ+ANNNI, and SC lattices, 
respectively. The renormalized coordination number for 
the SC lattice is significantly larger than the estimate for 
the Hciscnbcrg antiferromag nets, 3 C = 2.78. That is, the 
renormalized coordination number is different for models 
with different spin symmetry. 

The rapidness of the convergence of ((J 1 ) is quite re- 
markable. This makes a sharp contrast to the suscep- 
tibility of a single chain. For example, at J'/ J = 0.01 
the critical temperature, normalized by S 2 J, for S = 3/2 
ranges between 0.4 and 0.5. In this temperature range 
the single-chain susceptibility deviates from the asymp- 
totic scaling form by w 20% (see the inset of Fig- EJ) - On 
the other hand, the effective coordination number, which 
is merely the inverse of the single-chain susceptibility, di- 
vided by J', already converges well to the J' /J = limit 
(Fig. EJ . This implies that a cancellation of corrections 
to the asymptotic scaling behavior occurs at some deeper 
level. In other words, existence of a more fundamental 
equation than the modified CMF relation, i.e., Eq. (J2J 
with z replaced by £, is strongly suggested. 

Although the renormalization of the coordination num- 
ber is observed independent of S, the same as for the 
quantum Heisenberg antiferromagnets reported in Ref.y, 
one should however notice that £ has different values ev- 
idently for the SQ and SQ+ANNNI lattices, that is, it 
depends not only on the dimensionality but also on the 
real coordination number of the lattice, z. That means 
that the "universality" we observe in the renormalized 
coordination number is somewhat weaker than that for 
critical exponents, which depend only on the dimension- 
ality and not on the fine structure of the lattice, in the 
usual context of continuous phase transitions. In the next 
section we will unveil the origin of the "universality" as 
well as the "nonuniversality" in the renormalized coordi- 
nation number for the present Q1D Ising model. 



0.85 




0.4 1 1 — 1 — 1 — 1 

0.001 0.01 0.1 1 
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FIG. 5: J'-dependence of the effective coordination num- 
ber [Eq. |0] for the SQ, SQ+ANNNI, and SC lattices with 
S — 1/2 (squares), 1 (diamonds), and 3/2 (circles). The error 
bar of each point, except for the exact results for the S = 1/2 
SQ lattice, is much smaller than the symbol size. Horizontal 
dashed lines denote the limiting value for each lattice evalu- 
ated by taking an average of C(J') at J' = 0.001, 0.002, and 
0.005 for S = 1/2. 



in a transverse field and vice versa^* In this section, we 
first consider the mapping of the 5=1/2 classical Ising 
model. Generalization for S > 1 will be separately dis- 
cussed shortly. 

In order to derive an explicit relation between the pa- 
rameters of the classical and the quantum models, we 
again introduce the transfer matrix formulation for the 
d-dimensional classical Ising model. The partition func- 
tion of the S = 1/2 classical Ising model [Eq. JSJ)] is 
expressed as 

Z= lim tr{TD) M . (33) 

M— »oo 

Here T is a 2^= x 2^= matrix along chains and is a direct 
product of local transfer matrices 



V. MAPPING TO QUANTUM ISING MODEL 

A. Spin-1/2 case 

In the weak interchain coupling regime, J' / J -C 1, the 
correlation length along chains grows much faster than 
in the perpendicular directions as the temperature is de- 
creased, and it reaches the one-dimensional scaling re- 
gion, described by Eq. (|23|l , well before the correlation in 
the perpendicular directions starts to develop over more 
than a few times of the lattice spacing. In such a strongly 
anisotropic regime, it can be justified to take a continu- 
ous limit only in the chain direction, while keeping the 
lattice in the other directions discrete, without introduc- 
ing any uncontrollable systematic errors. In this limit, 
as is well-known, a d-dimensional classical Ising model is 
mapped onto a (d— l)-dimensional quantum Ising model 



e J/4T e -J/4T 
e -J/4T e J/4T 



(34) 



and D is a diagonal matrix representing the Boltzmann 
weight of each slice 



J' 



(35) 



On the other hand, the zero-temperature partition 
function of the quantum spin-i transverse- field Ising 
model defined on a (d — l)-dimensional lattice is writ- 
ten as 



Z a = lim tr e- H ^ T 



T^0 

lim 

A^0 



lim tre _AWt e _AWd e -AWt e _AWd 



(36) 
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where the quantum Hamiltonian 

7i q = Tit + TLd = — r — K o z r o z r 



(37) 



is defined on a (d - l)-dimensional lattice which has the 
same structure as a slice of the classical Ising model un- 
der consideration. Here er" (a = x, z) is a Pauli spin 
operator. 

By identifying T and D in the classical partition func- 
tion H33|) with e~ A7it and e _AWd of the transverse-field 
Ising model 136f) respectively, one can readily obtain the 
following relations: 



-J/2T 

AT 



tanh Ar and 
AK 



(38) 
(39) 



between the parameters of the two models, (J, J', T) and 
(K,T), or 



e -.//2T = tanh 



JT 



JT 



4ifT 4KT 



(40) 



by eliminating A in Eqs. (|38|l and (|39|l . The last equation 
is valid in the limit of T/J, J' /T — > 0, and relates the 
temperature T of the classical model in this limit to the 
strength of the transverse field T of the quantum model 
at the ground state. Especially, the critical points T c of 
the former at J' <C J and the quantum critical point T c 
of the latter is identified through 



-J/2T C 



AKT C 



(41) 



Now we are ready to calculate the effective coordina- 
tion number @ in the weak interchain coupling limit. 
From Eqs. © and l|4*T|) . 



C(J') 



1 



J' X ld (T c (J')) 
4r c JT C = £c 

J^o j' 4^ ~~ # 



42; 
J' 



-J/2T C 



(42) 



Thus the effective coordination number converges to a 
finite value in the J' = limit and its limiting value is 
exactly given by the quantum critical point, T c /K , of the 
corresponding (d— l)-dimensional quantum Ising model. 

One can indeed confirm that the analytic and numer- 
ical results for the renormalizcd coordination number, 
C = 1 and 3.0449(4) for the SQ and SC lattices, in the 
last section agrees completely with the quantum criti- 
cal point of the transverse-field Ising model, T c /K = 1 
(Ref.lH and 3.0440(7) (Ref-Hg) or 3.0450(2) (Ref.|13) on 
the single chain and the isotropic square lattice, respec- 
tively. Note that in the argument presented above, we do 
not assume any exact solutions except for the suscepti- 
bility of the genuinely one-dimensional model. Therefore 
Eq. (|42|l applies to arbitrary Q1D Ising models in any 
dimensions. (For example, it applies even to the models 
with frustration in the transverse directions.) 



B. General spin cases 

One may expect that a spin-S* classical Ising model 
could be similarly mapped onto a quantum Ising model 
with the identical spin size. However, this is not the case. 
Such a naive mapping, albeit formally possible, gives no 
physically meaningful weak interchain coupling limit. 

To see the difficulty for S > 1 more precisely, we con- 
sider the S — 1 case as an example. The 5 = 1 local 
transfer matrix along chains is given by 



( 



\ 



e J/T 
1 



-J/T 
1 



(43) 



-J/T l e J/T 



If this local transfer matrix is fit by an imaginary-time 
propagator exp(— AHt, r ) by the following S — 1 quan- 
tum Hamiltonian: 



(44) 



where the t" (a = x, z) is a spin-1 operator, T the 
transverse field, and D and E are the longitudinal and 
transverse crystal fields, respectively, we obtain 



Ar 



J 



-J/T 



V2T 
AD ~ J/T, and 
AE ~ -4J/T 



(45) 

(46) 
(47) 



up to the leading order in J/T. In the last equations, 
however, one should notice that Ar vanishes for T <C 
J (so that r itself can remain finite), whereas both D 
and \E\ diverge even after being multiplied by A. This 
desperate asymptotic behavior of the parameters of the 
mapped quantum system is due to the vanishing diagonal 
weight for the S z — state in the transfer matrix T r . 

Although for S > 1 the weight of the states with 
\S Z \ < S becomes exponentially small at low tempera- 
tures, disregarding such matrix elements in the transfer 
matrix will mislead us into an incorrect quantum spin 
Hamiltonian. In order to handle the anisotropic limit of 
the high-spin models with care, we introduce a unitary 
transformation defined below. 

As mentioned in Sec.[H] only the two largest eigenval- 
ues Ai and A2 of the local transfer matrix T r are domi- 
nant at low temperatures, and the other 2S — 1 eigenval- 
ues are exponentially small in comparison with Ai and 
A2. Now we introduce two column vectors v\ and V2, 
which are linear combinations of the eigenvectors u\ and 
1*2 associated with Ai and A2, respectively, 



= {u\ + u-i)/\Fi and 
V2 = (ui - u 2 )/V2. 



(48) 
(49) 



Clearly v\ and vi are orthogonal with each other by 
definition and furthermore they are interchanged by the 
spin inversion transformation, since u\ {ui) is symmetric 
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(antisymmetric) under the transformation. This implies 
that one can interpret V\ and t>2 as "wave functions" of a 
pseudo 5 = 1/2 quantum spin. Indeed, by a local unitary 
transformation by a matrix V r = (v%, ug, D3, . . . , V2S+1) 
with v\ and i>2 introduced above and Vi = Ui for i > 3, 
the local transfer matrix at low temperatures T -C J is 
reduced to a 2 x 2 matrix 



A+ a_ 

A_ A 1 



\ 



A2S+1/ 



(50) 



0/ 



where \± = (Ai ± A2)/2. The whole transfer matrix 
T = ® r T r of (25 + l) Ns x (25 + l) N ° is also reduced to 
a 2 s x 2 N " matrix by the global unitary transformation 
using V = ® r V r . We transform the diagonal matrix 
D at the same time to keep the partition function l|33(l 
invariant. Using Eqs. (|14|) - (|17|l . relevant matrix elements 
of D in this new basis are calculated as 

(V r ® V r+s )- 1 D r , r+s {V r <g> V r+S ) 

( e s\riT ^ 

e -S 2 J'/T 

e -S 2 j'/T 

P S 2 J'/T 







V 



(51) 



J 



The effective partition function of the spin-5 classical 
Ising model in the weak interchain coupling limit thus 
becomes equivalent to a spin-1/2 partition function, and 
then is mapped onto a spin-1/2 quantum Ising model, 
instead of spin-5, as the same as the 5=1/2 classical 
Ising model. Noticing that 



A+ A_ 
A_ A+ 



= exp 



exp 



%g(A + -A_) ilog£ 



A++A_ 



h l °gj£t log(A+-A_) y 



logAi ilogAi/A 2 
ilogAi/A 2 logAi 



(52) 



(53) 



one can explicitly write down the correspondence be- 
tween the spin-5 classical Ising model and the spin-1/2 
quantum transverse-field Ising model as 



log 



Ai(T) 
A 2 (T) 
S 2 J' 
T 



l/£i d (T) = 2Ar and 



= AK 



(54) 
(55) 



H, 3 
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FIG. 6: Scaling plot of the critical temperature for the SQ, 
SQ+ANNNI, and SC lattices with S = 1/2 (squares), 1 
(diamonds), and 3/2 (circles). Solid lines denote scaling 
function, 2ye~ 2 / y = Cx, with x = J'/SJ, y = T C /S 2 J, 
and C = limj^ C(J') [~ 1, 2.4755, and 3.0449 for SQ, 
SQ+ANNNI, and SC, respectively]. The error bar of each 
point is much smaller than the symbol size. The data for the 
SQ+ANNNI and SC lattices are shifted upwards by 1 and 2, 
respectively, for visibility. 



or 



V6d(T) = 



25VT 
KT 



(56) 



for J' <C T ^ J. Furthermore, using the low- 
temperature asymptotic form for the correlation length 
of a sing le chain [Eq. ffity]. we obtain a relation between 
the classical and the quantum critical point parameters 
as 



-2S 2 J/T C 



5JT C 
2T C K 



(57) 



which reduces to Eq. (|41|) for 5 = 1/2. The last equation 
suggests that two dimensionless variables, x = J' / SJ 
and y = T C /S 2 J, satisfy a scaling relation, 2ye~ 2/>v = Cx 
with an 5- independent constant C = T c /K . In Fig. 
this scaling relation is tested explicitly by using the 
Monte Carlo results in the last section. For each lat- 
tice the data in the weak interchain coupling regime 
J' /SJ < 0.1 collapse on a single curve as expected from 
the scaling form l)57p. 

Finally, using Eq. I|24|l . the effective coordination num- 
ber in the J' /J — limit is obtained as C = T c /K, 
irrespective of the spin size, which is a proof of the 5- 
independent renormalization of the coordination number 
suggested by the Monte Carlo simulation in the last sec- 
tion. One should note the nontrivial factor 1/5 in the 
scaling variable x ~ J' / SJ, which cancels with the same 
factor in the scaling form for the single-chain suscepti- 
bility [Eq. 124fl ] to make the effective coordination num- 
ber independent of 5. We also emphasize that our final 
result, ((J 1 ) — > T c /K, does not depend on the scaling 
conjecture for the susceptibility and correlation length 
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of a single chain; it can be deduced more directly from 
Eqs. |(iDJ and JSSJl instead. 

VI. SUMMARY AND DISCUSSION 

In this paper we presented our detailed analysis on 
the critical temperature and the effective coordination 
number of the Q1D general-spin Ising ferromagnets. The 
critical temperature is estimated down to J'/ J = 0.001 
with high accuracy by the cluster Monte Carlo method 
performed on infinite-length strips. It is demonstrated 
that as the interchain coupling is decreased the effective 
coordination number converges to a finite value which 
is definitely less than the coordination number z, quite 
rapidly. Furthermore, we show unambiguously that the 
renormalized coordination number is independent of the 
spin size S. This is the same as the observation for the 
Q1D and Q2D quantum Heisenberg antiferromagnets, re- 
ported in Ref. y, and thus is expected to be observed 
widely in other classes of quasi-low-dimensional magnets. 

We explain the origin of the renormalization of the 
effective coordination number by considering a mapping 
to quantum Ising models, where the spin-5 1 classical Ising 
model is rigorously shown to be mapped onto a spin-1/2 
quantum transverse-field Ising model, irrespective of the 
spin size S, in the weak interchain coupling limit, and 
the effective coordination number is given by the quan- 
tum critical point T c /K of the mapped quantum Ising 
model. It should be noticed that in the scaling relation 
of the critical temperature [Eq. (|57|) ]. a nontrivial fac- 
tor 1/5 appears, which cannot be deduced by a naive 
dimensional analysis. This factor cancels out with an- 
other nontrivial factor in the single-chain susceptibility 
[Eq. (23} ]) and then results in an S- independent scal- 
ing relation. Another remarkable fact to be noticed is 
that corrections to the asymptotic value in the effective 
coordination number are much smaller than those in the 
single-chain susceptibility itself, which suggests existence 
of a more fundamental theory, describing the critical tem- 
perature of strongly anisotropic magnets, than the mod- 
ified CMF relation, Eq. (5J combined with z = £. 

The renormalized coordination number is, however, 



found to be nonuniversal with respect to the lattice struc- 
ture. This is natural from the viewpoint of the mapped 
quantum Ising model, where the critical transverse field 
of a quantum critical point cannot be a universal quan- 
tity; it depends on the fine lattice structure as well as on 
the details of interaction. Thus an introduction of a weak 
further neighbor interaction or an anisotropy in slices, 
for example, will modify the renormalized coordination 
number. In this sense, the "universality," discussed in 
the present paper, is different from that asserted for crit- 
ical exponents in usual continuous phase transitions and 
thus inexact terminologically. 

Finally, we emphasize again that the weak interchain 
coupling limit of the Q1D classical Ising model is not 
the weak coupling limit ; it is still in the strongly corre- 
lated regime, i.e., in the vicinity of a quantum critical 
point, of the mapped quantum Ising model. This ex- 
plains why the CMF approximation remains inaccurate 
even in this limit. We expect a similar argument will ap- 
ply to the Q1D and Q2D quantum Heisenberg models. 
It is unclear and remains as an open problem, however, 
whether it is possible to generalize the present mapping 
approach to the quantum Heisenberg models, where a 
continuous axis, the imaginary-time axis, exists intrin- 
sically, and therefore extra care should be necessary to 
introduce another continuous axis representing a strong 
spatial anisotropy. 
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